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Abstract 

The mechanism of describing quantum states by standard probability (tomo- 
graphic one) instead of wave function or density matrix is elucidated. Quantum 
tomography is formulated in an abstract Hilbert space framework, by means of the 
identity decompositions in the Hilbert space of hermitian linear operators, with 
trace formula as scalar product of operators. Decompositions of identity are con- 
sidered with respect to over-complete families of projectors labeled by extra pa- 
rameters and containing a measure, depending on these parameters. It plays the 
role of a Gram-Schmidt orthonormalization kernel. When the measure is equal to 
one, the decomposition of identity coincides with a positive operator valued measure 
(POVM) decomposition. Examples of spin tomography, photon number tomography 
and symplectic tomography are reconsidered in this new framework. 
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1 Introduction 



Quantum states are described by vectors in a Hilbert space [1], or wave func- 
tions [2], in the case of pure states. In the case of mixed states, the density 
operators [3,4] are used instead to describe quantum states. On the other 
hand, different other tools have been introduced to describe quantum states 
by means of functions on phase space. These functions, like Wigner function 
[5], Husimi-Kano Q-quasidistribution [6,7], Sudarshan-Glauber diagonal (sin- 
gular) P-quasidistribution [8,9], contain informations about the quantum state 
which amount to the informations carried by the density matrix in an arbi- 
trary representation. In fact, these different quasidistributions are alternative, 
essentially equivalent forms of representing the density operators. These qua- 
sidistributions have some properties similar to those of classical probability 
distributions on phase space, but they are not fair joint probability distribu- 
tions since the uncertainty relation of position and momentum is incompatible 
with the existence of such probability density. 

Recently a tomographic approach to reconstruct Wigner functions from opti- 
cal tomograms was suggested [10,11]. The optical tomography approach was 
generalized to provide symplectic tomography [12,13]. In the tomographic ap- 
proach the quantum state is associated to a probability distribution depending 
on some extra parameters. This observation was used to develop a probabil- 
ity representation of quantum mechanics in which the tomographic probabil- 
ity distribution (tomogram) is considered as the primary object obeying an 
evolution equation of generalized Fokker-Planck type [14] and containing all 
informations on quantum state. Thus it is possible to formulate quantum me- 
chanics by describing a quantum state by fair probability distribution instead 
of wave functions or density matrices. By reading this chain of associations 
backwards, it is quite natural to ask if it is possible to provide an interpre- 
tation of tomograms directly at the level of the abstract Hilbert space. In 
particular this interpretation should work equally well for finite level systems 
(spin tomography) and generic systems. The interpretation can make obvious 
the mechanism of description of quantum states by fair probabilities instead 
of wave functions and density matrices. 

It is the aim of this paper to show how to provide such an interpretation with- 
out, however, indulging on more technical aspects (these will be considered 
elsewhere). The main idea consists of expressing the tomogram in terms of 
a scalar product in the space HI of rank-one projectors, that is in the linear 
space of operators acting on the space of quantum states. These projectors 
are connected with special families of vectors in the Hilbert space 7i of quan- 
tum states. The vectors are eigenvectors of families of operators depending 
on some extra parameters. We will obtain a decomposition of the identity op- 
erator in terms of a weighted sum of the projectors depending on the extra 
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parameter and determining the tomogram, so that any matrix and in partic- 
ular the density matrix can be obtained as a weighted linear combination of 
the basis vectors in the space HI. This explains why the inversion formula (re- 
construction formula) works for the tomographic maps. The tomograms can 
be constructed also for spin states [15,16]. The relation between tomographic 
maps and star-product quantization schemes was clarified in Ref. [17,18]. We 
develop our theory bearing in mind the case of spin tomography. However, 
that general picture of tomographic map is applicable to other kinds of to- 
mographies too, e.g. to photon number tomography [19,20,21] and symplectic 
tomography. To demonstrate the results we review the approach in which an 
n x n— matrix is considered as an n 2 — vector, used for instance in Ref. [22]. 

The paper is organized as follows. In the next section 2 we review the picture 
where a matrix is regarded as a vector. In section 3 we define the tomography 
in abstract finite dimensional Hilbert spaces and give new interpretation to 
the tomograms, both in terms of sets of rank-one projectors and of families of 
unitary, or Hermitian, operators. In section 4 we discuss in the light of our pic- 
ture some known examples as the spin tomography for the finite dimensional 
case, the photon number and the symplectic tomographies for the infinite di- 
mensional case, deriving new identity decompositions in terms of rank-one 
projectors in Hilbert space of bounded Hermitian operators. Conclusions and 
perspectives are discussed in section 5. 



2 Matrices as vectors 

In order to make clear how the tomographic approach provides relations con- 
necting probability distribution with density matrix elements we start with 
a very elementary example. We consider two Hilbert spaces Ti and HI. For 
simplicity we first identify the Hilbert spaces H with the qu-bit (i.e. spin 1/2) 
quantum state set, i.e. with vectors 

(1) 

Then the density matrix (density operator) for this pure state has the form 

It is well known (e.g., see [23]) that the set of operators acting on the Hilbert 
space H is a linear space. This space is a Hilbert space HI since one has the 
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scalar product of two operators A and B acting on the space Ti given by the 
formula 

(A\B) = Tr(A^B). (3) 

We use here the Dirac's notation for the scalar product. In fact, to write an 
operator (a matrix) as a vector, being very simple, is convenient. Thus the 
matrix of eq.(2) can be mapped onto a 4- vector using the rule 



Wi 



Ml 
Ml 
Ml 

Ml 



(4) 



This rule allows to reconstruct a matrix if the corresponding 4- vector is given. 
For example, given a 4- vector \A), one obtains the matrix A as: 



\A) 



ai 
a 2 
a 3 

04 



A 



O3 CZ4 



(5) 



Of course other "reconstructions" would be possible, and indeed this possi- 
bility has been exploited to consider alternative associative products on the 
vector space of matrices. [24] 

The scalar product of eq.(3) expressed in terms of matrices A and B is nothing 
but the standard vector scalar product given by 

{A\B) = (A\B) = Y t alb k (6) 
k=i 

The set of 4-vectors equipped with this scalar product is the Hilbert space H. 
Thus, having an initial Hilbert space of vectors TC of two-dimensional qu-bit, 
we have also the four-dimensional Hilbert space of 4-vectors EI = B(H), the 
linear operators acting on the Hilbert space Ti. The orthogonal basis in the 
space TL of spin up and down states (standard basis) 
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for instance, is complete. The completeness relation can be given in the form 
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of an equality valid in the Hilbert space HI = B(H), namely 

P 1 + P 2 = I H . 

Here the orthogonal projectors 



(8) 



Pi = |ei) (ei| 



1 




; Pi = |e 2 ) (e 2 | 




1 



(9) 



satisfying the relation (8) corresponds to the possibility to decompose each 
vector of Ti in terms of the basis vector \e\) and |e 2 ) , that is 



4>i |ei) + i> 2 |e 2 ) • 



(10) 



We considered in such details very simple properties of the linear spaces H and 
H, because essentially they are the basis of our considerations on the possibility 
to construct tomographic probabilities and to guarantee the existence of an 
inversion formula yielding the operator from its tomogram. In general, the 
completeness relation of a basis of H can be represented in the form 



The rank-one projectors 

h = m o*i ( 12 ) 

depend on a set of parameters ji (discrete or continuous, as well as finite or 
infinite) and, being in general non-orthogonal: 

P^Pjx' ± 0, 

they form a positive operator valued measure (POVM). 

On the other hand, in the four-dimensional space H one has the standard basis 
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, \B 2 ) = 
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\Bi) 

so that each vector \A) may be decomposed as 

\A) = "£a k \B k ). 

k=i 



(13) 



(14) 
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This decomposition of the vector \A) corresponds to the decomposition of the 
matrix A of eq.(5) in the form 





= ai 


1 


+ a 2 


1 


+ a 3 





+ a 4 





a 3 a 4 
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with respect to a basis of non-Hermitian operators. However, it is always 
possible to use a basis of Hermitian operators, for instance the unit matrix & 
together with the Pauli matrices a k) (k = 1, 2, 3). Then : 



1 
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(15) 



so that 

3 
k=0 

The relation between old and new components is 



(16) 







ai + a 4 


a 2 + a 3 




ai a 2 


1 


« 2 a 3 




i (a 2 - a 3 ) 


ai — a 4 




a 3 a 4 


" 2 



tto + «3 Oil ~ i a 2 



So, in general, if we were to use the coefficient of expansion in terms of a 
Hermitian set of n basis matrices, vectors corresponding to Hermitian opera- 
tors would be real and the Hermitian conjugation would act on them as the 
identity. In the case of non-Hermitian basis, those vectors are not real, in gen- 
eral, and belong to a real linear n 2 — dimensional submanifold of C n which is 
invariant under Hermitian conjugation. 

Now we transfer our discussion of the completeness relation of a basis of vec- 
tors, given for the two-dimensional Hilbert space Tl, to the case of the four- 
dimensional Hilbert space H. Thus we get 4-projectors 

F k = \B k )(B k \, fc = l,...,4. (17) 

The completeness relation in EI reads 

4 

E p * = fa, (18) 

k=l 
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where In is the four-dimensional unit matrix. Alternatively, in the four-dimensional 
space, one can have a POVM P^ with a set of parameters such that 

= (19) 

This property means that one can decompose any 2x2 matrix in terms of 
the 2x2 matrices corresponding to the projectors P^ . In matrix form the 
previous eq.(19) reads 

E„ ^»), h mn = SimSjn, i,j,m,n = 1,2. (20) 

That means that the index k in eq.(18) is considered a double index in labelling 
the elements of two by two matrices. In principle one can have more compli- 
cated conditions of completeness when the projectors P M are not orthogonal 
and have trace different from unity. Then a Gram-Schmidt orthogonalization 
procedure can be encoded by an extra kernel in the relation (20). We will see 
that precisely this situation takes place in some of the examples considered 
later on. 



3 The abstract Hilbert space definition of tomograms 

3.1 Tomographics sets 

In general, any kind of tomogram of a pure state is the positive real number 
W^(a,(3, ...), depending on a set of parameters (a, (3, ...) which label a set of 
states \a,P, ...) , defined as: 

W4a,f3,...) = \(a,[3,..W- (21) 

At a first glance, it seems quite difficult to read the tomogram as a scalar prod- 
uct, rather than a square modulus. Nevertheless, this is possible by thinking in 
terms of rank-one projectors rather than of (pure) states. In fact, by using the 
density operator p = and the projectors P a ,p,... = \a,(3,...) (a,(3,...\ , 

we may interpret the tomogram as a scalar product on the space of rank-one 
projectors: 

W^(a, [3, ...) = | (a, (3, ..#)| 2 = Tr {pP aA ..) . (22) 
This definition may be applied also in the case of an arbitrary density operator 
P- 

In the following we wish to characterize the sets of vectors \a,(3,...) which 
allow for a complete reconstruction of the state or an arbitrary density 
operator p, from the knowledge of its tomograms W^p(a, (3, ...). These sets will 
be called tomographic sets. In the light of our interpretation of the tomogram, 
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the meaning of such a reconstruction is nothing but a consequence of a decom- 
position of identity in the space of rank-one projectors in terms of the family 
\Pa,0,...) (Pa,p,...\ ? after taking into account that the projectors P a ,p,... m general 
are not orthogonal. We will discuss our interpretation in the finite dimensional 
Hilbert spaces. Our construction essentially goes along the following lines. 

Suppose a set {\e a /3)}a,i3=i °f n2 vectors of C n is found in such a way that 
the respective projectors \e a p) (e Q/ g| are a basis { |-Pfc) }^=i of C™ 2 = C n ®C n = 
B(C n ). We use here a collective index k instead of (cx,/3), e.g. k = (a — l)n+/3. 
By means of the Gram-Schmidt procedure, for instance, we may convert the 
basis {|P fc )}£ =1 into an orthonormal basis {|V})}™ =1 : 

n 2 

m = T,^\p k ) , {vm = k- (23) 
k=i 

In general, every element of the orthonormal basis {|V})} is a linear combi- 
nation of projectors, rather than a single projector like \Pk) associated to a 
vector of C n . 

There exists a decomposition of the identity on C n2 = B(C n ) 

9 9 9 

I n2 = £ |^.) {V 3 \ = £ 7 ; fc 7,7^Tr(A-) =J2K l Tr(P r ), (24) 

3=1 j,k,l=l 1=1 

where the Gram-Schmidt kernel Ki has been introduced 



Ki = E -ijClA (25) 
j,k=i 

We observe that Ki is a nonlinear function of the projectors Pk, because also 
the coefficients 7's depend on the projectors. 

We define the set {\e a f3)}a,f3=i °f n2 vectors of C™ a minimal tomographic set. 
The tomogram of a density matrix p with respect to this minimal tomographic 
set is defined by 

W p (a,[3)=Tr(\e a p)(e a p\p), (a, (3 = 1, n). (26) 

Then, from the decomposition of identity in terms of the tomographic projec- 
tors, we get an inversion formula for the density matrix p or any other operator 
on C n : 



P= E l* k l 3 iPiK{PkP) = E 

j,k,l=l V,v=l 



E ihiji{Pk)%Pi 
j,k,i=i 



(P)r- (27) 
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In other words, writing the previous equation in terms of matrix elements: 

(PV (28) 
we have the corresponding expression for the decomposition of the identity: 



(pW = E 

jX,V=l 



E ihiii{Pk)%{Pi)^> 
j,k,i=i 



n 2 n 2 n 2 

E ->J*7ii(A)^(^W = E(A)^(E i%i3i p i)^> 

j,k,l=l k=l 3,1=1 

= <V<W = (In <E> \)^,uu' ■ 



(29) 



A set containing more than n 2 vectors of C n is a tomographic set when it 
contains a minimal set. In other words, a tomographic set is such that any 
vector belongs to a (minimal) tomographic subset of n 2 vectors. In particular, 
a set is tomographic if any subset of n 2 vectors is a (minimal) tomographic 
set. 

Now, a basis of n 2 rank-one projectors can always be found. In fact, an or- 
thonormal basis of B(C n ) containing n 2 Hermitian operators is associated with 
the generators of the group U(n), multiplying each element by the imagi- 
nary unit %. Each generator, using its spectral decomposition, can be written 
in terms of projectors. Moreover, each projector can be expressed by means 
of rank-one projectors. So, from the spectral decompositions of the generators 
of U (n) , we may extract a basis of n 2 rank-one projectors. 

Alternatively, using an orthonormal basis {|a;)}™ =1 of n vectors of C n , we may 
define in B(C n ) an orthogonal basis of n 2 Hermitian operators given by: 



{(|a> (P\ + (a\) , i (\a) {(3\ - \(3) (a\)} , (a, (3 = 1, n). 



(30) 



Then, from their spectral decompositions we may extract a basis of n 2 rank- 
one projectors. 

For example, starting from a fiducial (orthonormal) basis of C 2 , two different 
suitable unitary operators 



—b* n* 



; a a a* a + b a b* a = 1; (a = 1,2) 



are needed to generate other two different basis of C 2 and obtain a tomographic 
set of six vectors containing three different minimal tomographic sets. In fact, 
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starting from the standard basis of C 2 , the matrix of 4- vectors 



1 a±al b\b\ a 2 a* 2 b 2 b\ 
— a\b\ a\b\ —a 2 b 2 a 2 b 2 
-a\b\ a\b\ -a*b* a* 2 b* 2 
1 b\b\ a\a[ b 2 b* 2 a 2 a 2 



has maximal rank, when the two different operators satisfy some extra condi- 
tion such as, for instance, 



^s(aibialbl) = det 



^ (cii&i) 3? (ai&i) 
3 (a 2 b 2 ) 3? (a 2 b 2 ) 



7^0. 



(31) 



This condition shows that the two complex numbers a±bi and a 2 b 2 cannot be 
proportional on the reals or, equivalently, they have different phases. 



Let us characterize more closely the manifold of rank-one projectors in the 

2 

real space M. n of Hermitian operators . Using the basis of the generators T k of 
U (n), with n — I and Trr fe = (k — 2, n 2 ) , we may express any Hermitian 
operator A as 

n 2 

A=5> fc r fc . (32) 

k=l 

The manifold of rank-one projectors is given by the vectors whose components 
{a k } G 1R™ 2 fulfill the conditions A 2 = A , TrA = 1 (which implies a 1 = l/n). 
In the dual space u* (n) of the Lie algebra of the generators r^, this manifold is 
an orbit O-p of the co- adjoint action of the group U(n). As the set of rank-one 
projectors may be identified with the projective space P (C n ) with 2 (n — 1) 
real dimensions, the orbit O-p is the only orbit with the same dimensions 
placed in the plane a 1 = l/n. Among the orbits, O-p is the one with lowest 
dimensionality, as the stability group of its points is just U(l)x U(n — 1). 
For instance, for n — 2, the rank-one projector set is the Bloch sphere S 2 : 
{ (a 2 ) 2 + (a 3 ) 2 + (a 4 ) 2 = 1/4} placed in the plane a 1 = 1/2. 

A minimal tomographic set therefore is a set of n 2 projectors {P(mk)}, where 
m k e Op, such that the linear span of {P{m k )} is the entire M™ 2 , that is 
det{K)'l + 0. Of course, Op is a maximal tomographic set. 
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3.2 Families of operators generating tomographic sets 



An interesting question is how to find a way to construct tomographic sets. 
This question may be answered in different but equivalent ways. We consider 
some of them here and provide a few well known examples to show how our 
proposal works. 

The first way consists in taking a fiducial rank-one projector P and acting on 
it with a suitable family of (at least n 2 ) unitary operators U a , depending on 
some parameters a. The family has to be chosen in such a manner that the 
set of projectors 



results into a tomographic one. This is granted only if the family U a is not con- 
tained in any proper subgroup of U (n) or, equivalently, if the group generated 
by the family U a is U(n). This condition is also sufficient if, moreover, from 
the family of unitary operators it is possible to extract, via the Cayley map, 
for instance, a basis for the Lie algebra u(n). Then, a family which is "skew" 
in the group U (n) is a suitable tomographic family of unitary operators. 

Alternatively, it is possible to start with a fiducial Hermitian operator A and 
to act on it with a "skew" family of unitary operators U a , generating a family 
of (iso-spectral) Hermitian operators 



Choosing Aq to be generic, i.e. with simple eigenvalues, the action of U a on 
the rank-one projectors associated with the eigenstates of A gives rise to a 
tomographic set of projectors. In other words, we may obtain a tomographic 
set from a suitable family of Hermitian operators. 

Taking U(n) as a tomographic family we obtain a (maximal) decomposition of 
the identity, analogous to eq. (24). However, integrating all projectors P over 
the symplectic orbit Op and using the volume f2 = oj n ~ l constructed with the 
canonical symplectic form u,we have to use a Hermitian kernel K{m), which 
is an operator valued function of the point m on the orbit O-p, that plays the 
same role of the Gram-Schmidt kernel in the minimal case: 



p a = u a p ui 



(33) 



A a = U a A Ul 



(34) 




(35) 



For instance, in the U (2) case, we have 




1 



l + 3cos# 3e _i *sin0 



(36) 
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so that, for any operator A, it results 

fllT f"K A 

A= / K{6,4>)Tt{P{6,(t>)A)smed6d(t>. (37) 
jo jo 

By using eq. (35) in the scalar product of any pair of operators A, B, we obtain 

Tr(AB)=Tr(A^ K(m)Tr(P(m)B)Q^ = Tr (^K(m)Aj Tr(P(m)B)Q 
= Tr (if (m)5) Tr(P(m)A)tt = Tr ^ jf P (m) Tr (i^T (m) 5) fi) . 

So, at least in a "weak sense", 

I n2 = / P(m)Tr(K(m)-)Q K(m) = [ P(m')Tr(K(m')K(m))n' . 

Finally, substituting the previous expression of K{m) in eq. (35), we obtain 
I n 2= f If Tr(K(m')K(m))P(m')n'] Tr(P(m)-)Q (38) 

JO v YJO v J 

in full analogy with eq. (24). Hence K(m) is just a Gram-Schmidt kernel, at 
least in a "weak sense". 

In general a tomographic family of unitary operators ranges between a minimal 
family and a maximal family, and is representative of the whole group U in) , 
so that necessarily the commutant of the family must be the identity. Then, 
a decomposition of unity must be available by integrating on the space of 
parameters of the family with a suitable kernel K{m). 



4 Examples 



4-1 Spin tomography 



For a qu-bit (i.e. a particle of spin 1/2), a tomographic iso-spectral two- 
parameters family of Hermitian operators is 



A(0,, 



cos 9 e ^sinO 
sin 9 — cos 9 



(39) 
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A given operator of the family corresponds to the component of the spin (up 
to a factor h/2) in the direction 

~n = (sin 9 cos 0, sin 9 sin 0, cos 9) . 

In fact 

It = A(0,<l>) (40) 
The spectrum of A(9, 0) and corresponding orthonormal eigenvectors are 



m± = ±1, \m + 9(f)} = 
while the respective projectors are 



e -^/2 cos #/ 2 



|m_#0) 



e -i^/2 s i n ^/2 
-e^/ 2 cos^/2 



(41) 



|m + #0) (m + #0| 



cos 2 |6> e"^ cos ±0 sin \9 



and 



sin ±0 cos \9 



sin 2 |0 



sin 2 \9 



-e _i * cos §0 sin \9 



\m-6<f>) {m 

They may be collected in the general form of a rank-one projector 



-e^ sin ±0 cos §0 



cos 2 \9 



P{0A) 



1 + cos 9 e _i * sin 9 
e lrt> sin 9 1 — cos 9 



(42) 



as 



P(9, 0) = |m+00) (m+00| ; P(vr - 0, tt + 0) = |m_00) (m_00| . 



Due to their form, two pairs of eigenvectors \m±9(f)) , \m±9'(f)') are not sufficient 
to yield a basis of projectors, so that at least three different operators of the 
family A(9, 0) are needed to construct a minimal tomographic set, as the 
previous example has shown. In the spin case, starting from the fiducial basis 
associated with A(0, 0) 



\m+) 



|m_) 



any A{9, 0)— basis is related to the fiducial basis via the unitary transformation 
f/(M): 

U{9^)\m ± ) = \m ± 94>) ) 
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where 



f/(M) 



-i<t>/2 



COS 9/2 e-^ /2 sinfl/2 



sinfl/2 



cos 0/2 



(43) 



Previous analysis has shown that two unitary operators U (9 , 0) , U \9' , 0') are 
a suitably "skew" set when ^s(e 1 ^ ~^ sin 9' sin 0) ^ 0. Thus, in the qu-bit case, 
a minimal set is obtained by taking a projector from each pair of eigenvectors 
and completing the basis with any other of the remaining projectors. So, as 
well known [27,28], only three independent directions of ~n are needed to 
reconstruct a spin 1/2 state, because the fourth projector is along any direction 
orthogonal to one of the first three directions. 

However, a decomposition of identity involving the whole family exists and is 
given by (see the previous eq.s (36), (37)) 



I 



2tt 



JO 



K(9, 0)Tr(P(#, 0)-) sin 9d9d<f>. 



(44) 



We observe that the kernel K(9, 0) is the simplest, the only one containing 
the same few spherical functions which appear in the projectors P(0,(f>), but 
it is not unique. In fact, a family of equivalent kernels can be obtained by 
adding to K(9, 0) any other kernel Ko{9, 0), containing only superpositions of 
spherical functions orthogonal to those of P(9, 0). 



As a matter of fact, in the qu-bit case the equation 

1= f 2 " [* P(8,<j))Tr(k(0,<l>)-)fin8ded<i> 
Jo Jo 



(45) 



holds in a strong sense. Hence K(9, 0) is just a Gram-Schmidt kernel in a 
strong sense. 

For the general case of spin —j<m< j, a decomposition of identity involving 
the whole family exists [18] and reads 

V / dtt (R(9, 0) \m) (m\ PJ(9, 0)) (K(m, 9, 0) W = 5(m'-s')5(m"-s") 

m=-j 

where R(9, 0) is a rotation through (9, 0) angles and dfl = sin 9d9d(p 1 while 
(X(m,0,0)W'=E E (2j 3 + l) 2 /(-l) m I>ffi(^fl,7)x 



(j 


3 h ^ 


' j j 


h | 




— m j 


K s' -s" 


m 3 J 



8tt 2 ' 



(46) 
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The problem of a minimal reconstruction formula for spin states was discussed 
in Refs. [27,28]. 



4-2 Photon number tomography 



This is an infinite dimensional case. However, the iso-spectral tomographic 
family of operators has a countable discrete spectrum n — 1, 2, 3, ...00, so that 
a generalization of our definitions is straightforward. We assume the fiducial 
basis {\n)} of the harmonic oscillator number operator a^a, and with the 
unitary family of displacement operators 

V (a) = exp (aa) — a*a*j (47) 

we generate the tomographic family A (a) depending on the complex param- 
eter a : 

A (a) = V (a) a ] aD ] (a) = (a f - a*) (a - a) , (48) 
and the respective basis of eigenvectors {|na)} = {V (a) \n)}. 

The photon number tomogram of a density operator p is 

W p (n, a) = Tr (\na) (na\ p) (49) 

while the inversion formula reads 

00 r rPn 

p=Y: —W p (n,a)KM(n, a ). (50) 

n=0 J 71 

The operator valued kernel is given by 

K^{n,a) = ^—(^-] n T{-a,-s), (51) 



l-s Vs-1 
where the operator T is 

T(a,s) = V(a)[^y a vHa). (52) 

Here s is a real parameter, — 1 < s < 1, which labels the family of equivalent 
kernels (n,a). 

The matrix form of eq. (50) in the position representation is 

d 2 a 



dx'dy' J (y' \na) (na\ x') (x\K^ (n, a) \y 



ln=0 



P (x', y') 
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To evaluate the matrix elements of the Gram-Schmidt kernel operator K^ s \ 
we first calculate those of the displacement operator T> (a). Remembering that 
a = (Q + iP) I y/2, we have 

V (a) = exp ((a — a*) -0= — i(a + a*) -j= j 



and putting 

we get 



a 



= {y- ip) /V2, 



(y\ exp (-i/iQ - ivP) \y) = S (y - y' - v) exp [i (-/ly' - fiv/2)] . 
Now the following map is useful 



s + 1 
s-1 



exp 



i In 



s + 1 
s-1 



where a determination of In has been chosen in such a way that t s > for 
s = 0. Then one readily obtains 



x\ exp 



-% (a f a + 1/2) T s + iT s /2~\ \y) 



= /n 1 ■ exp (i (x 2 + y 2 ) cot r s - + ^ \ 

^2m sm t s V L v 7 sin r s 2 J / 

and the matrix element (x\T(a,s) \y) results as 

/ dx'dy'(x\V(a) \x>) (x' \ ( i±±) \ y') (y'\V(-a) \y) = 1 , x 
J \s — 1/ \/27r«smr s 

exp f i ((x - z/) 2 + (y - ^) 2 ) cot r s - — — — - /i (x - y) + ^ J , 

\L smr s 2 J / 



so that eventually the matrix element ^1^^ (n,a) \yj reads 
x\K ( ' (n, a) \y 



1 — s \s — l) ^J2ixi sin r_ s 
exp ( % ((x + vf + (y + u) 2 ) cot r_ s - — + ^— ^ + ^ + // (x - y) + ^ ) 



(53) 



4-3 Symplectic tomography 



In the symplectic case (h = 1), we start from the fiducial basis of 
(improper) eigenvectors of the position operator Q : Q \X) = X \X) , whose 
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spectrum is the whole real axis: X e M. The two-real parameter family of 
unitary canonical operators S(p, u) : 

^(/i, z/) = expi-(QP + PQ)expi-(Q 2 + P 2 ); (/i = e x cos6,v = e" A sin6>) , 

generates both an iso-spectral family A(p, v) of Hermitian operators 

A(p, v) = S(p, u)QS\p, v) — pQ + vP 

and a tomographic set of (improper) eigenvectors \Xpv) = S(p, is) \X) , such 
that (X' pv\Xpv) = 5 (X — X') . In the position representation {\q)} it is, for 
v ^ 0: 



(q\Xpu) = (q\S(p,u)\X) = -j^=ex P 

V 2vr | | 



(54) 



Now, substituting the explicit expression of the symplectic tomogram : 

W P (X, /i, !/) = Tr (\Xpis) <*H p) = | (g |JS» </> (g'|/%> rfgrfg', 

in the well known inversion formula for (y\p\y'): 

p(y, y') = ^J v ) (y\ ex P (x - pQ- vp)\ W) dxdpdu, (55) 

we obtain 

x (</ |X/iz/) (X/iz/| g) dXdpdu} p(q, q)dqdq' . (56) 
In different terms 

1/'; Q, q) = J —dpdv (y\ exp [i (X - pQ - uP)} \y') (q f \Xpv) (Xpv\ q) 

= S(q-y)S(q'-y'), 

so that a partition of identity generated by the the tomographic set of rank-one 
projectors \Xpu) {Xpu\ appears in the inversion formula. 

Let us check the previous equation. From 

(y|exp [i (X-pQ-uP)] \y') 
= e lX (y\ exp (—iuP) exp (—ipQ) exp (—ipu/2) \y') 

= 6(y-y'-u) exp [i (X - py' - pu/2)} , (57) 
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we have the following expression of I(y, y'\ q, q') : 

dXdjjidu 



I(y, y'\ q, q') = J ^.X^ 5 (y - y - v) x 



(2tt)>| 



x exp [i (X - ny' - (iv/2)] exp i^-(q 2 - q' 2 ) - i^-{q - q') 



(58) 



Integrating over X we get \v\ 5 (q — q' — v) , which can be used to linearize the 
quadratic term: 



exp 



• A 4 / 2 /2\ 



5 (q — q' — u) = exp 



S(q-q'-v) (59) 



and we may write I(y, y'\ q, q') as 



i^(q + q'-2y'-u) 



(60) 



/ ~2n~ 6 ~y' ~ u ) 5 ^~ % ~ u ) ex p 

Integration over jj, yields 

I(y, y'; q, q') = J dvb (y - y' - v) 5 (q - q' - v) 25(q + q' - 2y' - v). (61) 
Eventually, we get the expected result: 



I(y, y'; q, q') = 25(q-q'-(y- y 1 )) S(q + q' - 2y' - (y - y')) 
= 26(q-y-(q'-y'))6(q-y + (q'-y')) 
= S(q-y)6(q'-y'). 



(62) 



We conclude this subsection recalling a problem posed by Pauli [25], wether it 
is possible to recover the state vector of a quantum system from the marginal 
probability distributions of the physical observables (e.g., position and momen- 
tum) of that system. The answer to the question is obviously negative [26]. 
The example of two squeezed states described by gaussian wave functions, 



ipi (x) = N exp —ax +i(3x , ip2 (x) = N exp — a*x + i(3x 



where > and (3* = (3, demonstrates readily this negative answer. The 
moduli of the functions are equal, and also the moduli of their Fourier trans- 
forms are equal. But the states are different since the scalar product of the 
two wave functions gives the fidelity which is not equal to one. So, different 
wave functions have the same marginal probability distributions of position 
and momentum. Now we are able to understand why the answer must be nega- 
tive: in fact, the family containing only the operators position and momentum 
is not tomographic because it is too small and cannot generate an inversion 
formula. 
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4-4 Squeeze tomography 



Finally, we discuss an example where an inversion formula is still lacking: the 
squeeze tomography [29] . The tomogram is defined using the same unitary op- 
erators S(fi, v) of the symplectic tomography, which acting on the fiducial basis 
{\n)} of the photon number tomography generate a basis of squeezed eigen- 
vectors {|n/xz/)} = {S(/i, v) \n)} of the squeezed tomographic family A sq (/i, v) 

A sq (fx, v) = S(fi, v)a)aS\ii, v) (63) 

In this case, however, the commutant of the family A sq (/i, v) contains the 
Parity operator and is nontrivial. Then the family is not a tomographic family 
strictu senso. Nevertheless, we get a true tomographic family by a restriction 
to the subspace of even wave functions. Then the existence of an inversion 
formula is granted. 



5 Conclusions 

We summarize the main results of the paper. The mechanism why quantum 
states can be described by fair tomographic probabilities instead of wave func- 
tions or density matrices was clarified. The mathematical reason for the pos- 
sibility to express the pure state projector \ip){ip\ (or density operator p) in 
terms of tomograms is based on the simple observation that any kind of to- 
mogram is just a scalar product of the projector, treated as a vector in the 
Hilbert space of operators, and a basis vector in this Hilbert space. The only 
property to be fulfilled is that the basis vectors in that Hilbert space are a 
complete (or even an overcomplete) set . For known examples of tomographies 
we have shown that it is always so. In view of this very elementary property it 
is even mysterious why the finding of the tomographic probability description 
of quantum states was done only relatively recently. 

Another result of the paper consists in finding explicit Gram-Schmidt orthog- 
onalizator kernels for symplectic and photon number tomographies. 

Consideration of tomograms of quantum states can be conceptually extended 
to the case of many degrees of freedom and even for the quantum field theory. 
In fact one needs only a pair of Hilbert spaces, H and B (H) , and constructing 
a tomographic basis in B (J~i). In principle, this consideration is extensible to 
infinite dimensional case (fields): one only has to add some extra ingredients 
to take into account the existence of different non-unique representations of 
the infinite Heisenberg-Weyl algebra and to use extra topological arguments. 
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Another obvious possibility consists in constructing the basis vectors to pro- 
vide tomographies by means of the eigenstates of quantum group operators. 
For the su q {2) case it only needs the introduction of operators dependent on 
the Casimir. 
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